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We investigate the time evolution of a generic and finite isolated quantum many-body system starting from 
a pure quantum state. We find the kinematical general canonical principle proposed by Popescu-Short- Winter 
for statistical mechanics can be built in a more solid ground by studying the thermalization, i.e. comparing 
the density matrices themselves rather than the measures of distances. In particular, this allows us to explic- 
itly identify that, from any instantaneous pure state after thermalization, the state of subsystem is like from a 
microcanonical ensemble or a generalized Gibbs ensemble, but neither a canonical nor a thermal ones due to 
finite-size effect. Our results are expected to bring the task of characterizing the state after thermalization to 
completion. In addition, thermalization of coupled systems with different temperatures corresponding to mixed 
initial states is studied. 
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I. INTRODUCTION 

Recently, studies of the foundation of statistical mechanics 
are enjoying a renaissance. Concerning about the canonical 
ensemble, it was pointed out that the equal a priori probability 
postulate of statistical mechanics [JJ, which is applied to the 
microcanonical ensemble constituted by all pure states satis- 
fying the energy constraint, is not necessary. Instead, it should 
be replaced bythe general canonical principle |01 or canoni- 
cal typicality yD. It states that the state of a subsystem can be 
obtained from a single pure state rather than a mixture of mi- 
crocanonical ensemble. Explicitly, suppose we have a weakly 
coupled system constituting by system 'S' and the heat bath 
'B', in case that the heat bath 'B' is large enough, starting 
from almost every pure state of the coupled system 'SB', the 
state of the system 'S' is approximately the canonical density 
matrix. The interest in these fundamental questions is also 
due to the fact that those postulates can be directly realized in 
recent experiments of ultracold quantum gases |4- 6]. Some 
related conclusions are obtained from different points of view 




The general canonical principle is of kinematic rather than 
dynamic [2^ ^ since it does not involve any time evolution 
of the state. Also the canonical density matrix from a pure 
state is because of the entanglement between system and the 
heat bath, otherwise the system will also be a pure state which 
cannot equal to a canonical ensemble. Dynamically, to result 
in a pure state between the coupled 'SB' system, the system 
as a whole should be an isolated system and also the initial 
state should be pure. And more importantly, all processes are 
of quantum mechanics such as the involvement of entangle- 
ment and Schrodinger equation for time evolution. Recently 
remarkable progress has been made in understanding the tem- 
poral evolution of the isolated many-body quantum system. 
In Ref.pi], Rigol et al. show that in a generic isolated system. 
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non-equilibrium dynamics is expected to result in thermaliza- 
tion: any pure states will relax to pure states in which the 
values of macroscopic quantities are stationary and universal. 
In particular, the thermalization can happen at the level of in- 
dividual eigenstates. However, the characteristic of the state 
after thermalization seems not clear, in particular about the 
relationship between the general canonical principle and the 
process of thermalization. 



In this paper, we will study a generic isolated quantum 
many-body system from the view points of both kinematics 
and dynamics. By using a generic thermalization model, we 
prove the general canonical principle numerically by showing 
that the state of 'S' from instantaneous pure state after ther- 
malization is similar as from a microcanonical ensemble. Our 
main method is based on comparing the exact form of den- 
sity matrices, rather than the measure of distance or the one- 
dimensional marginal momentum distribution, thus the results 
are more explicit and more precise. We then can clarify a cru- 
cial question that the state of a subsystem in dynamical equi- 
librium is like from a microcanonical ensemble or a general- 
ized Gibbs ensemble, but neither a canonical nor a thermal 
ones. Our conclusions are based on numerical results of a real 
physical system, so it is a numerical experiment to confirm the 
theoretical expectation. 



This paper is organized as follows: In Sect. II, we introduce 
the model. In Sect. Ill, we present definitions of several en- 
sembles of states. In Sect.IV, we obtain our main conclusions 
from numerical results. In Sect.V and Sect. VI, we show that 
our conclusions do not depend on specific lattice configura- 
tion and initial state. In Sect. VII, we study the thermaliza- 
tion of two coupled systems with different temperatures. In 
Sect. VIII we study the entanglement properties of the states 
in thermalization. Finally a summary of conclusions is pre- 
sented in Sect.IX. 
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II. MODEL 

Let us first present some fundamental concepts. Consider 
a generic quantum many-body system constituted by weakly 
coupled system 'S' and a large heat bath 'B', in thermody- 
namic equilibrium, the state of system is in a canonical form. 
In statistical mechanics to obtain this canonical state of 
'S', the composite system is described by the microcanoni- 
cal density matrix in which each pure state satisfying the con- 
straints of suitable total energy has equal probability. The state 
of the quantum system 'S', p^^\ can then be obtained by trac- 
ing out the freedoms of the heat bath from the microcanonical 
ensemble. By this process, the quantum state of the system 
in thermodynamic equilibrium, p''^', can be found to take the 
form of the canonical density matrix which is the thermal state 
defined as, 

nis) ^^-PH''^/z, (1) 

where /? is inverse temperature, H^^^ is the Hamiltonian of 
the system 'S', and the partition function is, Z — tie^^-^'^'^^ . 
The Hamiltonian of the system takes the form: 

H = -JJ2ib^k + h■c.) + uJ2f^^f'Jl (2) 

where indicates all pairs of nearest-neighbor sites on, a 
21 -site, two-dimensional lattice, see FiglTta) for its configu- 
ration, J is the hopping parameter, U is the nearest-neighbor 
repulsion parameter and we set it as, [/ = O.IJ. In Fi^TJ^)^ 
we consider a portion of the lattice, sites 1 to 8 (lower-right 
corner) as the system 'S', other parts of the lattice, sites 9 to 
21 (upper-left corner) are considered as the heat bath 'B'. In 
order to show that our results are general and do not depend on 
specific lattice configurations, we also extend the sites of heat 
bath 'B' from site 21 up to sites 25, see FiglHa). The coupling 
between 'S' and 'B' is through the interaction between sites 5 
and 12, iJ(-f), which can be turned on and off. So Hamiltonian 
^ can be rewritten as, 

H = H'^^'> + i? + i? . (3) 

As a whole, the coupled system 'SB' is an isolated quantum 
many-body system. In this paper, we consider that five hard- 
core bosons propagate in time on this lattice. 

The model ^ is in general non-integrable since the pecu- 
liar structure of the lattice in FiglTJa) has broken any spatial 
symmetry. By diagonalizing the Hamiltonian (|2]l with 21 sites 
(also up to 25 sites), we find that its energy level spacing dis- 
tribution is exactly a Wigner distribution [15.1, Pwigner{s) = 

^sexp(— ^^), see FiglTtb), which corresponds to noninte- 
grable system (or chaotic system). So thermalization is gen- 
erally expected to happen. 

m. TIME EVOLUTION AND ENSEMBLES 

We consider the initial state as |i/'(0)) which is the ground 
state of H'-^"> with all five hard-core bosons confined in eight 
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FIG. 1 : (a) Two-dimensional lattice on which five hard-core bosons 
propagate in time, the number of lattice sites can be selectively ex- 
tended from 21 sites as in Ref.|7] to 25 sites, (b) Level spacing dis- 
tribution of the energy spectrum of the Hamiltonian l|2j. The su- 
perimposed curves show a Poissonian (dashed green) and Wigner 
(solid red) distribution, characterizing integrable and non-integrable 
(chaotic) systems, respectively. 

sites of 'S', the interaction _ff is in off position initially. We 
can expand the initial-state wavefunction by the eigenstate ba- 
sis of the whole Hamiltonian H as, |'0(O)) = J2a C'qI^q), 
where Ca — {'^a\'4'{0)) are the overlaps between the initial 
pure state and the eigenstates |5'q), the corresponding eigen- 
values are Ea ■ 

Before proceed the dynamical relaxation, we consider dif- 
ferent sets of ensembles. The energy of the system which 
is conserved in the evolution can be found to be, Eq = 
{ip{0)\H\ip{0)) = — 5.06J. In statistical mechanics with the 
equal a priori probability postulate, the microcanonical en- 
semble can be defined as, 

Pmicro=-^ (4) 

where M is the number of Ea in the energy window, \Eo — 
Ea\ < Ai?, and AE — O.IJ is relatively small. The canoni- 
cal ensemble, which is the thermal state corresponding to H, 
is defined as, 

a 

with the probability distribution, P„ = e~^^°/Z, Z = 
J2a e"'^^°, /3 = l/ksT = 1/1.87 J, the inverse temperature 
/3 is consistently obtained from the equation, tr(pcano-ff) — 
e^^^" Ea/Z = En. We can also consider the general- 
ized Gibbs matrix jTlfTlll, 

Pgibbs=^|C„n«'a)(«'a|. (6) 

a 

We will mainly study the state of the quantum system 'S'. 
Corresponding to those ensembles, the density matrices of 'S' 
can be obtained by tracing out the freedoms of the heat bath 

(S) 

'B', p-^ — tiBRx, where X denotes the cases of micro- 
canonical, canonical and Gibbs, respectively. Together with 
thermal state fl^^^ defined in ([T]i, these four density matri- 
ces are schematically presented in Fig|2] From the statistical 
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mechanics, in thermodynamical limit with infinite large heat 



bath, as we mentioned. 



should take the form of thermal 



' ^micro ' 

State n'^^\ However, from Fig|2] it seems those two states are 

(S) 
■ 
micro 



not like each other. On the other hand, we may find that 

is more like Pgi^bs' there seems no other apparent similarities. 
Quantitatively, we can use a Hilbert-Schmidt norm to quantify 
the distance between two density operators, ||0|p = 
O is the difference of two matrices. We find. 



\\P, 



is) 



this confirms our observation from Figure|2] 



r micro "sibbs" 



(7) 



(8) 



Other distances are relatively larger ranging from 10 to 
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FIG. 2: Four density matrices of 'S' from different ensembles are 
presented. With Eq.([5]( as an example, the horizontal axes is Ea, the 
vertical axes is Pa in logarithmic scale, the results are almost linear, 
so slope corresponds to minus of inverse temperature. 

Now let's consider the dynamical relaxation. The time evo- 
lution begins with a sudden turning on the the interaction H^-^^ 
between sites 5 and 12. The initially confined five hard-core 
bosons begin to propagate on the whole lattice. The wave- 



function then evolves, \ip{t)) 



-iHt 



we have. 



(9) 



Our numerical experiment is to study the dynamical proper- 
ties of wavefunction in particular the time dependent 
reduced density operator of system 'S' compared with that of 
the states from various ensembles. In thermodynamic equilib- 
rium, we can also compare our result with that of the general 
canonical principle 101 • 



IV. THERMALIZATION AND GENERAL CANONICAL 
PRINCIPLE 

As mentioned in introduction, thermalization means that 
any pure states will relax to pure states in which the values 
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FIG. 3: (a) The time evolution of trace norm D{t) between p^-^\t) 
and thermal state The distances between thermal state with 

that of microcanonical and canonical ensembles are presented for 
comparison, (b) Time average of reduced density matrices 
for different scales of time period. All density matrices are almost 
the same. 



of macroscopic quantities are stationary and universal. By 
general canonical principle f?], it means that in thermody- 
namic equilibrium, from instantaneous pure state (wavefunc- 
tion) \'ij){t)), the state of 'S', p''^\t) equals approximately to 
thermal state 51*^'^) defined in ([T]). Numerically, we plot in 
Figl3ja) the time evolution of distance between two matri- 



ces. 



D{t) = \\p(s\t) 



For comparison, we also 

present the distances of p]^l^ro piano with fl^'^K As ex- 
pected, the time dependent distance quickly drops to almost 
zero. We confirm this result with different measures of dis- 
tance, the same conclusion can be made. Thus by thermaliza- 
tion, roughly, we confirm the general canonical principle by 
that any instantaneous pure state of 'SB' after thermalization 
leads approximately to the thermal state of 'S'. 

However, the following two important questions need be 
analyzed carefully by our numerical experiment: (i) Which 
state should be our goal state for comparing in time evolu- 
tion, thermal state O^'^^ or /ol^jcro^ (ii) What is the state of 
time average? Also we need to check what is the case for 

(S) 

different initial states. As we aheady found that p„jicro 
ferent from thermal state fl'-^^ in our system. Of course, it can 
be expected that those two states are the same if the freedoms 
of the heat bath 'B' are large enough. However for our finite 
two-dimensional lattice, the 13-site heat bath 'B' is not large 
enough compared with 8-site system 'S', so the first question 
arises. In the proof of general canonical principle in Ref.| J], 
both ensemble and time average is not necessary, only the dis- 
tance between reduced density operator of a typical pure state 
and the thermal state is used. Then in our system, what is the 
role of time average, thus the second question arises. We will 
clarify those questions by numerical results. 

From the time evolution of distance in Figl^a), we can 
guarantee that thermalization happens and the state is station- 
ary after t — 100. Obviously, we can find that the density 
matrices of time average for different scales of time period 
are almost the same, see Figl^b). Then it is reasonable to 
expect that after thermalization, the time average density ma- 
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FIG. 4: (a) The time average density matrix agrees with the instan- 
taneous density matrix, (b) The time average density matrix agrees 
well with the density matrix from generalized Gibbs ensemble. 



values. To show the results are independent of a wide range 
of initial states, we next will set the initial state differently, we 
will show that the conclusion still holds. Also we can consider 
different lattice configurations. 



V. THERMALIZATION WITH DIFFERENT LATTICE 
CONFIGURATIONS 

We expect that our results do not depend on specific lattice 
configurations. So we cast the heat bath 'B' with changeable 
number of freedoms, i.e. the sites 22-25 can be selectively 
added as part of the heat bath 'B', see Fi^a). Our main con- 
clusions expressed as a chain of equations always hold. 



trix is approximately equal to the instantaneous density matrix 
from a single pure state. 



it). 



(10) 



From the view point of distance, this can be confirmed in 
FiglSja). The comparison between time average density ma- 
trix with the instantaneous reduced density matrix is, to some 
extent, obvious and agrees well this result, see FigHJa). 

With the time evolution wavefunction (|9]l, the density ma- 
trix is simply, pit) = Ea,/3C':C;3e^(^=-^'')*|*„)(*^|. 
Considering thermalization happens which possesses a sta- 
tionary state, a long-time average of density matrix of 'SB' 
is expected to be. 



(11) 



This is the generalized Gibbs ensemble And the equiv- 
alence is confirmed by directly comparing those two density 
matrices, see Fig. |4|b). 

Remarkably, bear in mind that microcanonical ensemble is 
the same as the generalized Gibbs ensemble as shown in (l8|, 
we now reach a chain of equations from our numerical exper- 
iment step by step. 



AS) 



(12) 



Here we emphasize that this equation chain is not only based 
on the measure of distances, but rather based on the exact form 
of those density matrices which is, of course, more explicit 
and more precise. While for distance as in Fig|3ja), we can 
not even distinguish the cases of canonical ensemble and the 

microcanonical ensemble. 

I — t 

As in Ref.[2], next we should follow the standard statis- 
tical mechanics to obtain the thermal state from the state of 



microcanonical ensemble. The difference between p. 



(S) 



and 



thermal state fl^^^ should be a finite-size effect for our sys- 
tem. In this paper, we show that the equal a priori probability 
postulate can be replaced by a numerical experiment which 
confirms the general canonical principle. In addition, dynam- 
ically, our result is from a process of thermalization. We show 
any pure states relax to pure states corresponding to micro- 
canonical ensemble with stationary and universal macroscopic 



{S) 
^'gibbs 



(13) 



Fig-® is for case of 25-site lattice, all results are similar for 
cases with 21, 23, 24, 25 sites, respectively. 
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FIG. 5: Results with 25-site lattice. Four density matrices of 'S' from 
instantaneous pure state after thermalization, time average after ther- 
malization, the generalized Gibbs ensemble and the microcanonical 
ensemble are presented. They almost overlap with each other which 
confirms our conclusion. 

Next we present the states of 'S' for different lattice con- 
figurations, in particular, with number of lattice site of heat 
bath 'B' changing, see Figure Still thermal state and the 
case of microcanonical ensemble themselves are different, it 
seems that our system is not large enough. In principle, this 
difference is due to finite-size effect, in particular, our system 
is a general physical system with thermalization, the interac- 
tion between 'S' and 'B' is weak. So thermodynamical limit 
can be obtained when heat bath takes the limit to have infinite 
degrees of freedom. 

The state of 'S' from the instantaneous pure state is al- 
most stable and is approximately like from the the generalized 
Gibbs ensemble, the quantum fluctuation is almost negligible. 
One can find in Figure (|7J, states from instantaneous pure state 
in different time after thermalization are almost the same. 

So by those results, we find that our conclusions do not 
depend on specific lattice configurations, after thermalization, 
the state of 'S' is stable with negligible quantum fluctuation. 
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FIG. 6: Results for 21-site lattice and 25-site lattice. Three density 
matrices of 'S' are presented: time average after thermalization, the 
generalized Gibbs ensemble and the microcanonical ensemble. 
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FIG. 7: Results of 25-site lattice: density matrices of 'S' from in- 
stantaneous pure state at different time. 



VI. THERMALIZATION WITH DIFFERENT INITIAL 
STATES 

Our conclusion is that from any instantaneous pure state af- 
ter thermalization, the state of system 'S' is like from a micro- 
canonical ensemble. By definition, we know that the micro- 
canonical ensemble depends only on the energy of the initial 
state with a narrow energy window. So, we should expect that 
a wide range of initial states with fixed energy should relax to 
pure states with similar properties as the microcanonical en- 
semble. On the other hand, the state of time average is equiv- 



alent with the generalized Gibbs ensemble which depends on 
the distribution of Ca of eigenstates \'^a) for the initial state. 
It is thus initial state dependent. So it is necessary for numer- 
ical experiment to confirm our conclusion for different initial 
states. 

We consider the following six different initial states with 
approximately equal energy Eq. Still five hard-core bosons 
are on the lattice, however, initially the number of hard-core 
bosons in 'S' are: 5, 4, 3, 2, 1 and 0, and other hard-core 
bosons are confined in lattice sites of heat bath 'B', respec- 
tively. There is no interaction between 'S' and 'B' before 
time evolution, iJ'^^ , and we set the initial states as prod- 
uct (separable) pure states, {i/j'-^^O)) (g) |-0(-^)(O)). The en- 
ergies of 'S' and 'B' are, (^'(^^0)1^^^^ IV'^^HO)) = ^^o^^ 
and (?/;(^HO)|^^^^IV'^'^nO)> = E^o^\ respectively. How- 
ever, the total energy is {ip'^^^l {0)\ ® {ip'^^^0)\H\ip'^^^'^ (0)) (g) 
\^^^\0)) = E fa Eo, which is sUghtly different from the 
summation of energies of system 'S' and heat bath 'B'. This 
is because of the interaction H^-^\ which should be nonzero 
but must be weak in the thermalization. Explicitly, for our 
finite-size system, we have. 
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where S in table is the number of bosons initially in 'S'. We 
find the interaction between system 'S' and heat bath 'B' is 
weak, and so we have. 



Eq E ^ E, 



(S) 



E. 



(B) 



(14) 



The generalized Gibbs ensemble is defined by the distribu- 
tion of Ca which is the overlap between initial state and the 
eigenvectors \'^a) of the total Hamiltonian. Figure ([8]l shows 
the distribution of Ca for different initial states, so the gener- 
alized Gibbs ensembles for different initial states are almost 
the same, see Figure Q. 

To analyze our numerical data, we should confirm that our 
conclusion for different initial states should always hold. The 
numerical results are presented in Figure dTol ). here the time 
average is for short period of time so that it is approximately 
the instantaneous pure state. Since the microcanonical en- 
semble with fixed energy is independent of the specific initial 
state, the main point need be checked carefully is whether the 
case of microcanonical overlaps with the case of generalized 
Gibbs ensemble or the case of time average. We find that they 
all are almost the same. Thus we conclude that for a wide 
range of initially states, the state of 'S' is always like from a 
microcanonical ensemble. 



6 



°- 0.5 


^ 1 

^ 0.5 


„ 1 
LU 

^ 0.5h 


^ 1 
LU 

^ 0.5 




1 

0.5- 

1 

0.5- 



(a)B:0, S:5 


1 1 

1 — 


(b)B:1,S:4 


1 -=4^ — 


— r— \ \ 

\ (c)B:2,S:3 




— r— 1 1 

\ (d)B:3,S:2 ■ 


1 ==f— 

1 -4— 


— 1 — \ \ 

\ (e)B:4,S:1 ■ 

— ^ — \ \ 




(f)B:5,S:0 



-10 



a 10"' 



10"' 
10» 



a 10"- 



10"' 
10° 



a 10 ■ 



*:5» . • 



(a) B:0, S:5 



I, ^Ik I 

3 • 



(c) B:2, S 



(b) B:1,S:4 

-5 



5 



12 ' 



(d) B:3, S:2 



(e)B:4,S:1 ^Vt . ' 



(f) B:5, S:0 



FIG. 8: The distributions of Ca for different initial states, where 
initially the number of hard-core bosons on 'S' are 5,4,3,2,1 and 0, 
respectively 



FIG. 10: Three density matrices of 'S', which are time average after 
thermalization, the generalized Gibbs ensemble and the microcanon- 
ical ensemble, are presented. They agree well with each. 



FIG. 9: The generalized Gibbs ensembles are almost the same for 
different initial states, the number of hard-core bosons on 'S' are 
5,4,3,2,1 and 0, respectively. 
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FIG. 1 1 : With energies being almost the same, states from the mi- 
crocanonical ensembles overlap with each other. 



VII. THERMALIZATION OF TWO COUPLED SYSTEMS 
WITH DIFFERENT TEMPERATURES 



Figure (fTTI) shows that with energy fixed, the microcanoni- 
cal ensembles are almost the same, this is in accordance with 
the definition. 



In order to study the thermalization from different points 
of view, we next consider that the coupled systems 'S' and 
'B' initially are in different temperatures and the correspond- 
ing initial states are mixed states, this is in comparison that of 
previous sections, where the state of the whole system 'SB' 
is always pure. Since thermalization can happen in this cou- 
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pled system, a dynamical equilibrium can be reached and the 
whole system 'SB' should have a same temperature. Note 
that thermalization is closely related with quench where the 
Hamiltonian of the system changes quickly and the state of 
the system then evolves under the new Hamiltonian. 

In numerical calculations, we assume that, /3(^) = 1, 
= 2, which are inverse temperatures of systems 'S' and 
'B', respectively. We still restrict us to the case of five hard- 
core bosons, two of them are confined initially in 'S' and the 
remained three are in 'B' system. The initial state of 'S' is the 
thermal state as presented in ([TJ with (3 = 1, 

l]'(^Ue-^''\ (15) 

where the normalization of partition function is omitted with- 
out confusion. Similarly initial state of 'B' is also a thermal 
state, 

^,(B)^^-2H<-)^ (16) 

where /3'^' = 2. It is clear that both states are mixed. The 
thermalization begins with the turning on the interaction be- 
tween systems 'S' and 'B' by term H'-^K 

Figure (fTZt shows that the state of 'S' is initially a thermal 
state at time t = 0, then thermalization happens, the state 
changes at points t = 1,2,3. The dynamical equilibrium 
quickly reaches, we find that the state of 'S' at t = 7, 8, 9, 10 
is almost stable, which is shown in Figure (fT3l l. 
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FIG. 12: State of the system 'S' at time points i = 0, 1, 2, 3 starting 
from a thermal state. 

If we study the whole system 'SB', corresponding to differ- 
ent temperatures of 'S' and 'B' systems, the average temper- 
ature can be calculated as /3 = 1.4540. After thermalization, 
the state of 'S' is still different from a thermal state. Figure 
(fT4l i shows that the state of 'S' at time t = 10 is different from 
a thermal state though it starts from a thermal state. Here we 
argue that this is a finite-size effect, in case the heat bath 'B' 
is large enough, it will be like a thermal state. 

The involved states in this section are mixed states. Here we 
show that a dynamical equilibrium can be reached and ther- 
malization can happen under different conditions. 
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FIG. 13: State of the system 'S' at time points t = 7,8, 9, 10 start- 
ing from a thermal state. Thermalization happens and the instanta- 
neous states of 'S' are almost the same which shows the equilibrium 
reaches. 




FIG. 14: State of the system 'S' at time points t = 10 and the reduced 
density operator of 'S' (black) from a thermal state with /3 — 1.4540. 



VIII. THERMALIZATION AND DECOHERENCE 

In a generic isolated quantum system, thermalization of 
state 'S' which is initially a pure state, is induced by the in- 
teraction between the system 'S' and the heat bath 'B' which 
will destroy the coherence of the initial state. Thus entangle- 
ment in the initial pure state, in general, will decrease since 
of the decoherence from the heat bath. Still many aspects 
of entanglement are of interest in the time evolution in 'S' 
and 'B'. We use concurrence, one measure of entanglement 
|'22], to quantify the entanglement between different pairs of 
sites on the lattice, each is a qubit-qubit state. The results 
are shown in FiglTS] Here we would like to point out sev- 
eral interesting facts: (a) As expected, entanglement between 
nearest-neighbor sites, C4,7, is relatively larger, while in time 
evolution, it does not decrease monotonically which indicates 
a non-Markovian dynamics, (b) With thermalization, all pair- 
type entanglement converge to almost zero including the sites 
in heat bath 'B' and the coupling states of (5,12), though the 
whole system is always a 21 -partite pure state. This may indi- 



8 




0.4 






0.2 


(c) 


^13.14 












' ■ ^13.17 










t 



FIG. 15: The time evolution of entanglement measured by concur- 
rence C for pairs of sites. For example, C4.6 is tlie concurrence 
between sites 4 and 6. 



cate a property of thermal state which can be described classi- 
cally so without entanglement, (c) The entanglement of pairs 
in 'B', is relatively large for a period in the beginning of ther- 
malization while decrease finally to zero. This indicates the 
quantum aspects of the thermalization for an isolated quantum 
system, it also indicates that entanglement may be necessary 
in the process to relax to a stationary state. 

Besides pair-wise entanglement, it is also of interest to 
study the local von Neumann entropy as measure of entan- 
glement between coupled systems. Since initially the state 
of 'S' is a pure state, when the interaction between 'S' and 
'B' is turned on, the von Neumann entropy of 'S' increases 
quickly and reaches almost a constant. This is also an indi- 
cation of thermalization, i.e. thermalization can be measured 
by chaotic of system 'S' quantified by local von Neumann en- 
tropy. We remark that since the overall state is pure, this von 
Neumann entropy is a measure of entanglement between sys- 
tems 'S' and 'B'. 



IX. CONCLUSIONS AND DISCUSSIONS 

In conclusion, by numerical experiment, we prove a chain 
of equations and reach a conclusion that after thermalization, 
the instantaneous density operator of 'S' is equivalent to the 
state from a microcanonical ensemble. Thus dynamically, we 
prove the general canonical principle which can replace the 
equal a priori probability postulate for microcanonical ensem- 
ble in statistical mechanics. We use the density matrix itself 
rather than distance to analyze our results thus their deriva- 
tions are more explicit and precise, we identify that the dif- 
ference between state from the microcanonical ensemble and 
the thermal state is a finite-size effect. In the generic isolated 
quantum system, we systematically present various ensembles 
and clarify their relations numerically. In addition, we show 
in the thermalization, though the state of 'SB' is always a pure 
state, its pair-type entanglement of system 'S' is almost zero 
like a classical state, but it jumps from zero in the beginning 
of thermalization to finite value and drops again to zero when 
thermal equilibrium is finally reached. This shows that entan- 
glement may be useful for thermalization. 

The main method in this paper is to compare density matri- 
ces schematically, however, the result agrees well with method 
of calculating the distances evaluated by such as Hilbert- 
Schmidt norm as used in Ref.|2]. It is already shown in 
Ref.|7] that the one-dimension momentum distributions are 
close to that of microcanonical ensemble. Since expectation 
values of observables are based on measurement operators ap- 
plied on density matrix, the method by comparing directly 
density matrices seems more complete and explicit. 

The main conclusion of this paper is based on a real phys- 
ical system in which thermalization happens. The process of 
thermalization is realized by a numerical experiment so that 
the result does not depend on any assumptions. Considering 
that the system is generic, our conclusion should be general. 

We would like to thank useful discussions with C. P. Sun 
and H. Dong. This work is supported by NSFC and "973" 
program (2010CB922904, 201 1CB921500). 
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